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Splitting Instability of a Multiply Charged Vortex in a Bose-Einstein Condensate 
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We consider the splitting mechanism of a multiply charged vortex into singly charged vortices in 
a Bose-Einstein condensate confined in a harmonic potential at zero temperature. The Bogoliubov 
equations support unstable modes with complex eigenfrequencies (CE modes), which cause the 
splitting instability without the influence of thermal atoms. The investigation of the excitation 
spectra shows that the negative-energy (NE) mode plays an important role in the appearance of 
the CE modes. The configuration of vortices in splitting is determined by the angular momentum 
of the associated NE mode. This structure has also been confirmed by the numerical simulation of 
the time-dependent Gross-Pitaevskii equation. 

PACS numbers: 03.75.Kk, 05.30.Jp, 67.40.Vs, 67.40.Db 



I. INTRODUCTION 

The discovery of Bose-Einstein condensates (BECs) of 
alkali- metal atomic gases enables us to make fundamental 
investigations into superfluidity and quantized vortices. 
Vortices in BECs have been created in three ways: phase 
imprinting 0, 0] , optical spoon stirring Q and topologi- 
cal phase engineering 0, [5j . By using the third method, 
multiply charged singular vortices are created in BECs, 
which have not been achieved in other systems such as 
superfluid 3 He, 4 He 0,0. 

Recently, MIT group succeeded in imprinting vor- 
tices in Bose-Einstein condensates by using topological 
phases They observed doubly and quadruply quan- 
tized vortices in a non-rotating trap. Although such a 
multiply charged vortex is expected to be unstable and 
split into singly charged vortices, no obvious splitting has 
been observed in the experiment. Unfortunately, since 
the center of the external potential changes from a min- 
imum to a saddle in the vortex imprinting process, the 
vortex state cannot be trapped. Therefore, the fate of 
the vortices has not been revealed in the experiment. 

In a non-rotating trap, the energy of a vortex state in- 
creases in proportion to the square of the winding num- 
ber [8|. When there are some dissipative processes such 
as scattering with thermal atoms, the multiply charged 
vortex should split into singly quantized vortices and es- 
cape from the condensate. The existence of the negative 
energy eigenvalues among the excitation spectra shows 
that the vortex state continuously turns into a vortex- 
free state without any energy barriers 0, IToL ITU . Then, 
what happens in the system in which all atoms are con- 
densed and scattering with thermal atoms can be negli- 
gible? Does the vortex remain in the condensate without 
splitting? 

There is another candidate for splitting instability. 
For a multiply charged vortex in a BEC, the Bogoli- 
ubov equations possess complex eigenfrequencies in cer- 



tain regions of the parameter space of the interaction 
strength [l2l Il3l Il4| . This fact means that there is 
an eigenmode growing exponentially and the conden- 
sate becomes unstable against infinitesimal perturba- 
tions. Therefore the vortex state is intrinsically unstable. 

In this paper, we consider the origin of the complex- 
eigenfrequency (CE) modes and their effect upon the 
vortex splitting process. A vortex with a winding num- 
ber L > 2 may split into various states, while a doubly 
charged vortex merely splits into a pair of singly charged 
vortices 0]. We will treat a quadruply quantized vor- 
tex in concrete calculations in this paper. In Sec. [HI 
we briefly introduce the formalisms based on the mean 
field theory. The orthonormal conditions for CE modes 
are constructed there. In Sec. IIIII a clear insight into 
the origin of CE and negative-energy (NE) modes is pro- 
vided by the collective excitation spectra of the Bogoli- 
ubov equation. The possible structures of splitting are 
also considered in Sec. lIVI showing the results of numer- 
ical simulations of the Gross-Pitaevskii equation (GPE). 



II. FORMALISM 

A. Stationary state with a multiply charged vortex 

BECs of weakly interacting atoms are well described 
with the GPE given by 



(1) 
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where ^(r, t) is a condensate wave function, M is the 
mass of the trapped atoms, /i is the chemical poten- 
tial, and g = Anh 2 a/M represents the strength of the 
interparticle interaction. We treat repulsively interact- 
ing atoms and assume the s-wave scattering length a 
to be positive. An external trap potential has an axi- 
ally symmetrical form with trap frequencies u r and lo z 
as Vtr(r) = \Mu? r {x 2 + y 2 ) + \Muj 2 z z 2 . For simplicity, 
we assume a pancake-shaped BEC, i.e. uj r -C ui z , and 
fix the z-dependence of the orderparameter as \&(r, t) = 



2 



tp(x,y,t)ip z (z), where ip z (z) is the normalized ground 
state wave function of a one-dimensional harmonic oscil- 
lator. Then we treat the two-dimensional wave function 
tj)(x,y,t) which is normalized to the total atom number 
N. 

It is convenient to introduce the dimensionless time 
and space variables t = uj r t and (x,y) = (x,y)/anc>i an d 
the normalized order parameter tp — anQt^/y/N, where 
ano is the harmonic oscillator length: oho = \J fr/2Muj r . 
Energies are also scaled by huj r , for example, jl — fj,/Hu> r . 
Then the GPE is rewritten in a simple form 



(2) 



where V 2 = d 2 + <9 2 , r 2 = x 2 +y 2 ,rj = 8waN and we omit 
tilde for simplicity. 

Assuming an axial symmetry around a vortex line, we 
express the condensate wave function with an L-charged 
vortex in equilibrium as 

Mr,0) =A(ry Le , 1 = 0,1,2,..., (3) 

where A(r) is a real function and satisfies the equations 

[C(L)+ V A(r) 2 ]A(r)=0, (4) 

(5) 



. . d 2 1 d L 2 1 , 
^L)^-—--- + -r + -r 2 



dr 2 



dr 



The chemical potential /i is found from the normaliza- 
tion condition 2it rdr {A(r)} 2 = 1. We assume L to 
be positive without loss of generality. The results of nu- 
merical calculations for L = 4 are presented in Sees. IIIII 
andllVl 



B. Collective excitations 

To study the collective excitations of a BEC, we add 
small fluctuations to the stationary state as ip(r,6,t) — 
[A(r) + f(r, 9, £)] e tLe ', where / is small and complex. As- 
suming that the excitations are periodic in 9 with period 
27r, we expand / into a Fourier series: 

f(r,9,t)= ]T [ Ul (r,ty w +v:(r,t)e- m ]. (6) 

1=0,1,2,- 

Substituting ip(r, 9, t) to Eq. J2J) and linearizing with re- 
spect to /, one obtains the well-known Bogoliubov equa- 
tion 



H;w/(r) = wjawi(r), (7) 



where 



-*>-($>)■ - ( o - 1 ; ■ < 8 » 

and Ti.i is a symmetrical matrix 

# _ f C(L + 1) + 2t 1 A 2 rjA 2 , 

ni ~\ r,A 2 £(L-l) + 2r 1 A 2 > ' (Jl 



We are interested in the eigenstates of afii and take 
the time dependence of the excitations as w; (r, t) — 
wj(r) exp(— iuiit). It should be noted here that the Bo- 
goliubov equations with different I are independent. 

Usually, the summation in Eq. JBJ is done for both 
positive and negative I which correspond to the angu- 
lar momenta of excitations, and the normalization con- 
dition (|u;| 2 — \vi\ 2 ) = 1(> 0) is imposed, where (• • •) = 
2ir J Q ■ ■ ■ rdr. In this case, however, we remove this nor- 
malization condition and restrict I > 0, since the eigen- 
value equation for — / is quite the same as those of I 
when one replaces ui,vi,u>i with and — o>_;, re- 

spectively. Then we define the angular momentum l cx of 
excitations by comparing the amplitudes of components 
exp(±iZ0) as 



Zex = I X Sgn(|li;| 2 - H 2 ). 



(10) 



-I, 



To specify the sign of the angular momenta, we use w 

and wj v ' for the eigenfunctions of l eK = I and l cx = - 
respectively. The excitation energy is also defined by 

e, =uj 1 x (M 2 - \ Vl \ 2 ), (11) 

reflecting the relation uji = —lo-i- This is confirmed by 
expanding the Hamiltonian up to the second order in 
perturbations. 



C. Complex eigenmodes and orthonormality 

Equation (JJJ is allowed to have complex eigenfrequen- 
cies since it is a non-Hermite eigenvalue equation |l2l . 
Here, we derive the normalization and orthogonality con- 
ditions for eigenmodes including CE modes. 

Consider two eigenmodes w/ m .„ of Eq. J7J with eigen- 
frequencies u>i m . n . Using the property of the Hermitian 
operator Hf. (wf n Tliwi m ) = (wf m HiWi n ) , where T de- 
notes the transpose, one obtains the equation 

(ujlm - wj„)(wf„ffwi m ) = 0. (12) 

When two eigenstates have different eigenfrequencies 
(u>l m ^ win), they are orthogonal: (w^<7W; m ) = 0. Espe- 
cially, in case when w;^ is complex with a complex eigen- 
value (dip, its complex conjugate w* is also an eigenstate 
with eigenvalue lu^ . The orthogonality condition for w;^ 
and w ; *^ is written as 



l 2 )=0. 



(13) 



According to Eqs. ljTU|) and lfTT|) . a CE mode corresponds 
to the excitation with zero angular momentum and zero 
excitation energy. 

Equation (|12|l indicates that (w^crw;„) can be nor- 
malized. For a real-eigenfrequency (RE) mode, we choose 
the phase so that the eigenfunction is real. Then the nor- 
malization constant is defined to be consistent with its 
angular momentum as 



/ (u)T~ (u), 
( W L °" w ln) 



1, 



, (v)T« (v) v 



1. 



(14) 
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For a CE mode, we divide the eigenfunction into two real 
functions w^ R ''' as w^, = -l=(w( R ) + rw^). The phase 

of wj m is determined in such a way that the following 
conditions are satisfied, 



(w< R > T Jw( R ) 
(w (R)T- w (IK 



1, 



1. 



(15a) 
(15b) 



With these conditions, the CE mode indeed satisfies 
Eq. I|13[) and the normalization conditions (w^erw/p) = 

(i' w fu) T v^lu) = 1- According to Eq. Ijl5a|l . the func- 
tion w' R ''^ correspond to positive- and negative-angular- 
momentum excitations (not eigenstates) . We use the ex- 
pression w|^, andwj^, instead ofw( R ''), respectively, and 
treat them and other eigenmodes equally. 

In this way, the orthonormal basis set w ; v " 
constructed with the following conditions: 



(u,v) 



has been 



>) T ; 
f lm 



(u)\ 
In I 



= 5n 
= 0. 



/ (v)T- (v)> 



(16a) 
(16b) 



III. ENERGY SPECTRUM 

A. Collective excitations with r\ — 

First, we investigate at non- interacting limit [r] = 0). 
The equilibrium state satisfies C(L)A(r) = 0. The eigen 
function F% and eigenvalue E% are given analytically by 



Eh = |Z|+2n + l, 



-r 2 /4 



(17) 
(18) 



where n = 0, 1, 2, • • • denotes the radial quantum number 
and &n( x ) 1S * ne generalized Laguerre polynomial func- 
tion. The coefficient C% is defined by the normalization 
condition (F^{r) 2 ) — 1. The stationary state Ao(r) and 
its chemical potential /zq are 



A)(r)=F L (r), mo 



L + l. 



(19) 
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FIG. 1: Excitation spectra of a quadruply charged vortex 
state at non- interacting limit. The symbols o and x represent 
and £;„ , respectively. 



Eq. fTTJl , the excitation energy of is £, 



-ft 



while that of W ; ( ^ is = il^' . Figure [T] shows the 
excitation spectra with rj = with respect to the equilib- 
rium state with a quadruply quantized vortex (L = 4). 
There are a few negative energies among the spectra of 
1 < I < 2L — 1. These NE modes come from the fact that 
-00 is not a ground state. For a non-interacting BEC, 
NE modes correspond to the eigenstates in the trap with 
lower energies than ip . The existence of these various 
NE modes shows how unfavorable the multiply charged 
vortex is. 

In the case of a non-isolated system, the energy of the 
condensate does not need to be conserved. It decreases 
through the growth of NE modes, leading to the vortex 
decay. In an isolated system, however, the energy must 
be conserved and the excitation which can cause the vor- 
tex decay is not a NE mode but a CE mode, as we will 
show in the following. 



,(u) 



,(u) 



(v) 
In ' 



The Bogoliubov equation with rj = decouples into 
the independent equations 



C{L- 



m 



-n t v h 



(20a) 
(20b) 



where we write eigenfunctions and eigenfrequencies with 
rj = in capital letters. The solutions are Ui = 
F r ^ +l ,Vi = F%~ 1 , i.e. the set of eigenmodes is given by 



W 



(u) 



F, 



L+l, 



r) 



W 



(v) 



(21) 



B. 77-dependence of energy levels 



To calculate the energy spectrum for 77 > 0, we expand 
the eigenmodes w/ in Wj"' v ^ as 



N 

71 = 



(22) 



having eigenfrequencies ft 



ft 



M 



(u) 



\L + l\ + 2n - L and 



\L — l\ — 2n + L, respectively. According to 



Substituting Eq. to Eq. 0, we rewrite the 

Bogoliubov equation as eigenvalue equations for 



340 730 T\ 

FIG. 2: Imaginary parts of eigenfrequencies for (a) I = 5, 6 
and (b) / = 2, 3, 4 as functions of the interaction strength r\. 
The result for I = 6 is magnified a hundred times. The CE 
mode at P differs from others with I = 4 as to the associated 
NE mode (see text). 

(ao,Po, ■ ■ ■ , «jv,/3jv): 

+(E% +1 - ii - ui)a n = 0, (23a) 

N 

r, £ [{F£- l A 2 F£- l )a m + 2{F^ 1 A 2 F^~ l ) f} m ] 

m=0 

+(E%- l - t jt + wi)0 n = O, (23b) 

where n = 0, 1, • • • , N and A(r), are the numerical so- 
lutions of Eq. (g}. 

We have numerically solved Eq. (|23|l for L = 4 and 
< I < 8 in the range of < r] < 4000. The number 
of terms summed in Eq. (|23|l is N = 50. The NE modes 
exist for 1 < I < 7 as expected from the result of rj = 0. 
The CE modes appear for 2 < I < 6. The imaginary 
parts of complex eigenfrequencies are shown in Fig. [21 
which agrees well with the results of doubly and triply 
charged vortices [l2^ . 

The origin of the CE mode is made clear by plotting 
them as functions of rj. Figure |3| shows the excitation 
spectra with 1 = 2. Each eigenstate is identified by a 
radial quantum number n = 0, 1, 2, • • ■ and a sign of an- 
gular momentum u,v. The solid lines in Figs.|«Ja) and 
(b) represent e^n an d e 2 n> respectively. There is a NE 
mode wf, \ and — 62^(> 0) is plotted with a broken line 
in both figures. We also plot the real parts of the complex 
frequencies with dotted lines. It is clear from Fig.[3]that 
a conjugate pair of CE modes appears instead of two RE 

modes and in the regions where e^„ + 4ia = 0- 
As we mentioned in Sec. [H] a CE mode corresponds to 
an excitation with zero excitation energy and zero an- 
gular momentum, which are made up of two excitations 




FIG. 3: Excitation spectra for I = 2. The solid lines cor- 
respond to excitation energies with (a) positive angular mo- 
menta and (b) negative angular momenta. In both figures, 
the negative energy eigenvalues and the real part of the com- 
plex eigenfrequency are plotted by broken and dotted lines, 
respectively 

such as -^=(w^ ± iv^2o)- The real and imaginary parts, 
which are divided so as to satisfy the orthonormality, 
have a physical meaning: they correspond to positive- 
and negative-angular-momentum excitations, which con- 
tinuously change to RE modes as rj changes. 

The second quantized description simplifies this pic- 
ture [T2I H*5| . A boson field operator is defined as 
ijj(r) = ipo(r) + 4>(r), where the c-number function 
tpo ( r ) denotes the condensate wave function and (j)(r) 
is the fluctuation part. We decompose cf>(r) as cf>(r) = 
e iLe Ex,iJu%\ry ie axi n Un\r)e- a6 a{j. Here, a[ ln 

is a creation operator associated with a fluctuation wj^ . 
Let us fix the interaction strength, for example, to rj — 
200 where a pair of the CE modes -^(yv^o ^* w 20 ) ex ists. 
Then the Bogoliubov Hamiltonian is written as 

A=u,v (i,n)^(2,0) 

+ Reuj c (a\ 20 a u2 o - aJ 2 o a v2oJ 

+ Imw c (al 20 al 20 + a u2 oa V 20 ) , (24) 
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1000 2000 3000 4000 



FIG. 4: Excitation spectra for / = 4. The description of lines 
are in Fig. |3Ja) . There are two NE modes and both turn into 
CE excitations. The point P corresponds to P in Fig. [5] (sec 
text). 

where uj c is the complex eigenfrequency, and we have used 
the relations 

Rec c = (w^ T H 2 w^) = -(w< v )T 7i 2 wM), (25a) 
Ima; c = (wg T H 2 wM). (25b) 

Equation (|24|l clearly shows that the two excitations 

^ are coupled to each other and should be created or 
annihilated together. 

In the above discussion, the NE mode is important for 
the appearance of CE modes. In the case of I = 2, there 
exist only one NE mode, • Therefore, all CE modes 
with I = 2 are found by tracing e^o • We also present the 
spectra for I = 4 in Fig. 0] where the solid lines represent 
£4^' . There exist two negative energy eigenvalues e^j and 
e 41 ^ which are plotted by broken lines. We have confirmed 
that the NE mode , as well as w^q , turns into CE 

modes in the region where + = is satisfied (at 
point P in Figs. [5f b) and^J. We have obtained all of the 
CE modes by inspecting the behaviors of NE modes (see 
Fig.©. 

The ^-dependence of negative energy eigenvalues is 
also important in this mechanism. Since their depen- 
dence differs from that of positive energies, many CE 
modes appear (Figs. EJa) and^J. The qualitative ac- 
count for the ^-dependence can be given by introducing 
the effective potential 

V eS (h r) = -r 2 + + 2 V A(r) 2 - fi, (26) 

4 r z 



which is easily derived from the Bogoliubov equation for 
negative-angular-momentum (— I) excitations. When the 
equilibrium state contains a vortex, the A 2 term produces 
a potential well centered at the vortex core ^(J- O n the 
other hand, the second term represents the centrifugal 
potential depending upon the angular momentum of the 
fluctuation. This term makes the well narrower as \L — l\ 
increases. A NE mode corresponds to a bound state of 
the well, and therefore shows different 77-dependence from 
others. 

There are two effects on the negative energy eigenval- 
ues associated with the change of 77. By increasing 77, (i) 
the chemical potential /i increases as well, leading to a 
deeper well and smaller energy eigenvalues and (ii) the 
coherence length (i.e. vortex core size) decreases, lead- 
ing to a narrower well and larger energy eigenvalues. In 
the case of I — 4(= L), the negative energy eigenvalues 
decrease rapidly as 77 increases by (i). In contrast, / = 1 
negative energy eigenvalue increases due to (ii). There- 
fore, the condition + ej$ = cannot be satisfied with 
any eigenmodes w^, and no CE mode appears. 

IV. PATTERNS OF VORTEX SPLITTING 

So far our discussion has been based on the spectrum 
analysis. We have found that the CE mode corresponds 
to a zero-energy excitation and should grow in an isolated 
system. In this section we consider the possible patterns 
of vortex splitting. 

A. Vortex structure with elementary excitations 

Here, we analyze the structure of BEC with fluctua- 
tions. The wave function perturbed by an eigenmode 
W/ is written as ip = [A(r) + 5 (uie iW + vfe~ m )] e iLe , 
where 8 is a small and real constant. 

The structure of the vortex core in equilibrium can be 
described with an asymptotic form A(r) ~ r L (r — » 0). 
The eigenmode at the center also behaves as ui ~ rl i+/ l 
and 17 ~ r \ L ~ l \ _ Then the perturbed wave function is 
written in the lowest order of r as ip ~ S rl i_ 'le^ L_ ^ e , 
i.e. an (L — /)-charged vortex exists at the center. One 
also finds that the condensate has /-fold symmetry by 
linearizing the density with respect to 5: 

\ip(r, 9)\ 2 = A 2 + 2AS [Re(u/ + vi)cos(W) 

-Im(ui + vi)sm(W)]. (27) 

Assuming the conservation of the total angular momen- 
tum, splitting of an L-charged vortex obeys the rules: 
(i) an (L — /)-charged vortex stays at the trap center, 
and (ii) I singly quantized vortices are arranged in i-fold 
symmetry. 

Although the above rules are applicable to both RE 
and CE modes, only the latter can grow exponentially 
and cause large changes of the condensate with a small 
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perturbation. In the following, we investigate the effect 
of CE modes to the time development of a vortex state. 



B. Projection to elementary excitations 

We solve the GPE numerically and simulate the time 
development of splitting process. Regarding the time 
evolution of a condensate as individual developments 
of elementary excitations, we again expand the time- 
depending wave function as 



0.02 



<Kr,0,*) = lo(t)A( r y Le 
+ 51 MM)e 



w 



vUr,t)e~ lW ] e lLe .(28) 



Conversely, ui and vi are written by using ?p(r,9,t) as 



ui(r,t) = 
vi(r,t) = 



fvM,i) e - 4(i+0e , (29a) 

Z7T 



<m 

2^ 



r(r,9,t)e 



i{L-l)6 



(29b) 



We expand w; (r, t) in eigenmodes of Bogoliubov equation 
as 



w;(r, t) 



E 

A=u,vra=0,l,2, 



E 



(A) 



(*)w& } (r). 



(30) 



By using the orthonormal conditions l|16|) , the coefficient 

is given by ^' v \t) = ±(w[^ v) (r)<7w ; (r,i)), where the 
upper (lower) sign is for a superscript u (v). 

In the extent of the linear approximation, 



Tin 
>,v) 



for 



a.s 



a RE mode oscillates with the eigenfrequency uj\ n 
7/n' V ' ) x ex P( — * w ;n' V ^)- Then the amplitude is constant. 
For a conjugate pair of CE modes -j= (w£$ ± «w|^, ) with 
eigenfrequencies uj c and w*, their time dependence are 
given by ^(7^ ± Him') a exp[-i(Re u c )t =p (Imw c )f], 
one of which is growing and the another is diminishing. 



C. Time development of vortex structures 

The solution of Eq. Q remains stationary without any 
perturbations. We distort the trap potential in the l v - 
fold symmetrical form as V tI (r, 0) — \r 2 [1 + <5cos(Z v #)] 
during a short period. By choosing the symmetry of the 
distortion, the eigenmodes of I = l v are selectively ex- 
cited. 

The distortion is small enough not to cause the distin- 
guishable splitting of the multiply charged vortex. When 
all of the excited modes have real eigenfrequencies, the 
configuration of cores does not change in time. In our 
simulation, no obvious splitting occurs at 77 = 500 where 
no CE mode exists. 

We have investigated how CE modes affect the config- 
uration of vortices. We have simulated at r\ — 730 where 



:o.oi 



1 1 

A (distortion off) 
\ 


(a) 


■] plus mode 






minus mode 



10 15 
t [trap unit] 



20 




100 150 
t [trap unit] 



200 



FIG. 5: Time development of elementary excitations, (a) 
Initial changes. The broken lines correspond to the amplitude 
of the CE modes (plus and minux modes), which coincide with 
the solid lines corresponding to the functions exp(±Imu; c t). 
The amplitude of other modes plotted with the dotted lines 
are constant, (b) Further time development. The plus and 
minus modes (the solid and broken lines, respectively) interact 
with each other and oscillate. The amplitude of other modes 
(dotted lines) are also affected by the oscillation and increase 
more and more. 



exists a pair of CE modes -^(wgg ±iw^ ), and distorted 
with l v = 3. Figure f5] shows the time dependence of 
73 n , which is the typical behavior of excitations includ- 
ing CE modes. During all of the amplitudes of eigen- 
modes are small, their time-dependence agrees well with 
the linear approximation as shown in Fig.rSJa): the am- 
plitude of one of the CE modes (plus mode) increases as 
exp(Imw c i), that of another CE mode (minus mode) de- 
creases as exp(— lmw c i), and those of RE modes are con- 
stant. Further time development is shown in Fig. f5Jb). 
When the plus mode increases beyond the linear approxi- 
mation, it stops growing and the minus mode starts grow- 
ing instead. They interact with each other and begin to 
oscillate. Other excitations are also enlarged more and 
more throuth a non-linear process. Here, RE modes of 
I = l v are selectively enlarged. (In the case of / v = 2, we 
found the RE modes of I = 4 are also excited through 
this non-linear process.) 

Figures Ufa), (b) and (c) are the images of conden- 
sates at the points A, B and C in Fig. |3] respectively. 
The initial distortion causes little change of the distribu- 
tion (Fig.HJIa)). As the complex mode grows, the vortex 
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(a) (b) (c) 




FIG. 6: Contour plots of density profiles at the points A, B 
and C in Fig. (a) A multiply charged vortex just after the 
distortion, (b), (c) Vortices arranged in three- fold symmetry. 
Vortex cores move outward (b) and into the center (c) in the 
oscillation. 



split into four singly quantized vortices arranged in three 
fold symmetry (Fig. EJb)). The interaction of the CE 
modes results in the oscillation of distances among the 
cores. The divided cores in Fig. EJb) move toward the 
center again as shown in Fig. EJc) and oscillate between 
these two states. In the oscillation, singly charged vor- 
tices never return to the quadruply charged vortex. Their 
average distance becomes larger and larger. Simultane- 
ously, the center core begins to vibrate and condensate 
becomes distorted more and more. The energy of the 
multiply charged vortex is gradually transfered to the 
excitation energies of RE modes. 

We have also investigated the coexistence region of sev- 
eral CE modes. We have simulated at r\ = 340 where CE 
modes of I = 2,3 and 4 are exist. Figures|7fa) , (b) and (c) 
are the images of splitting after distorted with l v = 2, 3 
and 4, respectively. Vortices in each image are arranged 
in Z v -fold symmetry. This fact means that the splitting 
patterns can be controlled by the initial distortion. 

As we have mentioned before, the symmetry of the 
configuration of vortices agrees with Z v , i.e. the angular 
momentum of excited states. At the center of the trap, 
however, a doubly quantized vortex, which appears in 
I = 2 splitting, seems unstable and soon decays into two 
singly charged vortices (Fig. Ufa)). 

V. CONCLUSION 

We have calculated the collective excitation spectra of 
a BEC with a quadruply charged vortex, and discussed 
the origin of CE modes. A CE mode is decomposed into 
two excitations, whose total energy and angular momen- 
tum are equal to zero. This fact allows the CE mode to 
grow without any dissipative processes. To satisfy the 
energy conservation, one of the elements must be a NE 



mode. Its ^-dependence affects the appearance of the CE 
mode. We have found all of the CE modes existing in the 
region < T) < 4000 (Fig. |£J|. 

The possible patterns of vortex splitting are classified 
by the angular momenta of fluctuations. These struc- 
tures have been confirmed by the numerical simulations 

of splitting. Moreover, the selection of the patterns has 
(a) (b) (c) 




FIG. 7: Images of splitting patterns at r\ — 340 after trap 
distortion with (a) l v = 2, (b) l v = 3, and (c) l v = 4. The 
uppers are density profiles and the bottoms are corresponding 
phase profiles. Three CE modes are selectively excited by 
choosing l v . 

found to be possible by the manipulation of the trap 
symmetry. The intrinsic splitting are caused by a CE 
mode, and if one choose the r\ (i.e. the number of trapped 
atoms) so that a CE mode exits, the splitting correspond- 
ing to the CE mode occurs. The CE mode is efficiently 
excited by controlling the symmetry of the trap, even in 
case when several CE modes coexist. 

In this paper, we have assumed a pancake-shaped BEC 
and taken a 2-dimensional approach. In a cigarette- 
shaped BEC as in the experiments, 3D simulation shows 
that the splitting cannot be recognized by time-of-flight 
because of the excitation along a vortex line Q . 
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